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We show that a gravitational interaction between the derivative of the Ricci scalar curvature and the baryon- 
number current dynamically breaks CPT in an expanding universe and, combined with baryon-number- 
violating interactions, can drive the universe towards an equilibrium baryon asymmetry that is observationally 
acceptable. 
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The successful predictions of Big-Bang Nucleosynthesis 
(BBN) |1], highly precise measurements of the cosmic mi- 
crowave background |2], and the absence of intense radia- 
tion from matter-antimatter annihilation [3] all indicate that 
the universe contains an excess of matter over antimatter. Nu- 
merically, the baryon-to-entropy ratio is ns/s — 9.2^q^ x 
10 -11 . What remains a mystery is how the baryon asymmetry 
was generated. 

The contemporary view is that the baryon asymmetry is 
generated dynamically as the universe expands and cools. 
Sakharov |4j] argued that three conditions are necessary: (1) 
baryon-number non-conserving interactions; (2) C and CP 
violation; and, (3) a departure from thermal equilibrium. To 
satisfy the latter two conditions, the conventional approach 
has been to introduce interactions that violate C and CP in 
vacuo and a period when the universe is out of thermal equi- 
librium. 

In this paper, we propose a mechanism that generates an ob- 
servationally acceptable baryon asymmetry while maintaining 
thermal equilibrium. The key ingredient is a CP-violating in- 
teraction between the derivative of the Ricci scalar curvature 
TZ and the baryon number (B) current J^: 
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where M* is the cutoff scale of the effective theory. It is not 
necessary that J M be the B current; any current that leads to 
net B — L charge in equilibrium (L is the lepton number), so 
that the asymmetry will not be wiped out by the electroweak 
anomaly H, is sufficient for our purpose. It is natural to ex- 
pect such an operator in the low-energy effective field theory 
of quantum gravity if M* is of order the reduced Planck scale 
M P = (8-kGn)- 1 / 2 ~ 2.4 x 10 18 GeV. We also note that it 
can be obtained in supergravity theories from a higher dimen- 
sional operator in the Kahler potential. 

The interaction in Eq. Q violates CP and is CPT conserv- 
ing in vacuo. However, this interaction dynamically breaks 
CPT in an expanding universe and biases the energetics in 
favor of causing an asymmetry between particles and anti- 
particles. 

To generate a baryon asymmetry using the interaction in 
0, we also require that there be B-violating processes in 
thermal equilibrium. We denote the temperature at which B- 
violation decouples by Tq. Given these ingredients, the B 
asymmetry in our setup is generated as follows. In an ex- 
panding universe, where TZ ~ H 2 and TZ are non-zero (where 



a dot means time derivative), the interaction in Eq. Q gives 
opposite sign energy contributions that differ for particle ver- 
sus antiparticle, and thereby dynamically violates CPT. This 
modifies thermal equilibrium distributions in a similar fashion 
as a chemical potential p ~ ±TZ/M 2 , driving the universe to- 
wards non-zero equilibrium B asymmetry via the B-violating 
interactions. Once the temperature drops below Trj, as the 
universe expands and cools, a net asymmetry remains: 
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Our approach is closely related to "spontaneous baryoge- 
nesis" [6], which relies on the derivative coupling between a 
spatially uniform scalar field and the baryon number current, 
(d^tpjJ^. With tp, though, the construction is considerably 
more complicated. The scalar, essential to this mechanism, 
has to be added by hand, whereas the term in Eq. is ex- 
pected to be present in an effective theory of gravity. The 
initial conditions for p must be separately specified and jus- 
tified: ip must be forced to evolve homogeneously in one di- 
rection versus the other to produce an asymmetry and must be 
spatially uniform. In contrast, the time-evolution of TZ oc H 2 
is required in a cosmological background and it is highly spa- 
tially uniform because the universe is highly homogeneous. 
The oscillation of ip around its minimum is also a compli- 
cation because the average (p is zero, tending to cancel the 
asymmetry |7], whereas the mean value of TZ ~ H 3 does not 
vanish. 

To produce baryon asymmetry by the gravitational interac- 
tion in Eq. Q, several factors have to be considered. For a 
constant equation of state w, where w is the ratio of the pres- 
sure p to the energy density p, the scalar curvature is propor- 
tional to (1 — 3w), and its time derivative is given by 
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We will examine B generation for some cosmologically im- 
portant values of w. 

We start from the radiation-dominated era following in- 
flation when w w 1/3. If w were equal to precisely 1/3, 
then the right-hand side of Eq. Q would vanish and there 
would be no effect. However, w = 1/3 only applies in the 
limit of exact conformal invariance, T^ — 0. In practice, 
interactions among massless particles lead to running cou- 
pling constants, and, hence, the trace anomaly that makes 
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T£ oc f3(g)F^F f _ l „ ^ 0. The thermodynamic potential of 
a plasma of the SU(N C ) gauge theory, with coupling g and Nf 
flavors, has been worked out in detail 1 8]. It leads to 



l-3w = 
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up to 0(g 5 ) corrections, where the last factor in the numer- 
ator is the beta function coefficient. Typical gauge groups 
and matter content at very high energies can easily yield 
1 — 3w ~ 10 _2 -10~ 1 . There may also be mass thresholds 
that lead to conformal violation. Then, Eq. Q gives 
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The upper bound on tensor mode fluctuations constrains the 
inflationary scale to be Mj < 3.3 x 10 16 GeV 0, and 
obviously for this scenario Tp, < Trd < Mi, where T R p> 
is the temperature at which the universe becomes radiation 
dominated (i.e., the reheat temperature in this case). It is re- 
markable that the asymmetry can be sufficiently large even for 
M* ~ Mp if Td — Mi. In this case, the scenario predicts 
that tensor mode fluctuations should soon be observed. 

A second case of cosmic relevance is w = 0. This cor- 
responds to the matter domination epoch which character- 
izes, for example, conventional perturbative reheating via a 
scalar field <j> osc , as it oscillates around the minimum of a 
quadratic potential. During the oscillation phase, p oc a~ 3 , 
a oc i 2 / 3 , and <j) osc decays at a rate V into radiation, whose 
energy density becomes equal to that of the scalar field when 
H ~ T ~ T% D /M P . Therefore, 



This result is 3-4- orders of magnitude enhanced relative to 
Eq. and allows for T RD ~ 10 14 GeV. 

A third possibility is to generate the baryon asymmetry 
while a non-thermal component with w > 1 /3 dominates the 
universe. The non-thermal energy component decreases more 
rapidly than radiation, so there is no need for this component 
to decay into additional radiation and produce more entropy 
in order to enter the radiation-dominated epoch. Hence, once 
ns/s is set at Tp, it remains constant. The absence of fur- 
ther dilution opens up the range of allowed parameters. One 
example where this can occur is in the ekpyrotic |9] or cyclic 
ficll universe where the kinetic energy density of a scalar field 
<f> dominates immediately after the bang over a smaller, sub- 
dominant radiation component. The scalar field is the mod- 
ulus that describes the interbrane separation. Another pos- 
sibility is that the inflaton <f> falls down a steep potential at 
the end of inflation and shoots out as a massless scalar field 
flUl . More general w can be realized, for example, by a co- 
herent oscillation of <f> about the minimum of the potential 
V{4>) = \4> 2N /Mp N ~ A with a coupling constant A. This form 
of V((f>) yields w = (N - 1)/(N + 1), where 1/3 < w < 1 
for N > 2 fl2ll . Such a potential is natural in supersym- 
metry using a discrete symmetry, and it is easy to verify that 
soft supersymmetry breaking effects do not spoil its desired 
behavior. 

We begin our analysis with the universe dominated by 
p<t> ~ a~ 3l - 1+w \ which decreases faster than a sub-dominant 
radiation component p R ~ a~ 4 . The ^-dominated universe 
expands as a oc t 2 /[ 3 ( 1+tu )], while the temperature drops as 

T(t) = T RD [a RD /a(t)] = T RD (t RD /t) 2 /^ 1+w ^ . There- 
fore, 



Pos 



PR 



T, 



II D 



1 RD 



a R D 
a 

a RD \ 3 / 2 



(6) 
(7) 



1 RD 



a R D 



3(l+w) 



_ „ 4 fa RD \ 
p R ~ T RD (— j 



(11) 
(12) 



The latter equation suggests T ~ T R o(a R D/a) 3 ^ 8 - At the 
time of decoupling Td > T R d and Eq. (|2} gives 
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This asymmetry, however, is diluted by a continuous produc- 
tion of entropy. The dilution factor is given by (T R d/Td) 5 
and, hence, the final asymmetry is 
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Within the linear approximation made in Eq. (|2), the initial 
asymmetry in Eq. (|8j cannot be larger than 0(1). Therefore 
T R d cannot be smaller than about 10 _2 Td to obtain the cor- 
rect baryon asymmetry, which gives an upper limit 
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Combining these relations, we find the asymmetry 
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Since we have the freedom in this scenario to make Td sig- 
nificantly greater than T R d and since w > 1/3, the baryon 
asymmetry in Eq. Jl 31 can be significantly enhanced relative 
to that in Eq. © by a factor {Td/Trd) 3 ^- 1 ^ 2 . Henceforth 
we focus on the case w > 1/3. 

Next, we discuss the origin of the i?-violating interaction 
that is necessary for any of the baryogenesis scenarios con- 
sidered here. To keep the discussion general, we assume 
that £?-violating interactions are generated by an operator 
Ob of mass dimension D = 4 + n. The rate of such in- 
teractions is given by T R = T 2n+1 /M B n , where M R is 
the mass scale associated with Ob- Decoupling of B vio- 
lating processes occurs at T ~ Td, when T R falls below 
H - (T 2 D /Mp)(T/T RD )^ 1+w ^ 2 . The decoupling temper- 
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ature Tjj is then estimated to be 

When w = 1 and n = 1, Eq. (I14> is not applicable since the 
interactions with dimension five operators, such as those re- 
sponsible for the neutrino mass seesaw mechanism are 
either in thermal equilibrium all the time or never. Using 
Eqs. O and <03}, we can identify the values of (Trd, Mb) 
that generate ng/s ~ 1CP 10 . Assuming n = 1, Fig^shows 
the result for various values of w. The plots for n > 1 are 
similar. 

In supergravity theories, gravitino production places severe 
bounds on the highest temperature T max attained in the early 
universe. Apart from order unity coefficients, the gravitino 
abundance Y3/2 is given by 

^- n -f^°-'w(M 1 '"""- <i5> 

where we have used Boltzmann's equation dn 3 / 2 /dt + 
3Hn 3 / 2 = a eff n R' w i tn (Jeff ~ 1/Mj, and n# ~ T 3 . 
We set T maa; = Tjj for numerical estimates below. The 
bounds come from two constraints: (?) ensuring that the prod- 
ucts of BBN will not be dissociated by late gravitino decays 
and (z'z) avoiding overclosure of the universe by gravitinos, if 
they are the Lightest Supersymmetric Particles (LSP's), or by 
the LSP's that would be produced in gravitino decays. Here, 
we assume that the gravitino has a mass TO3/2 ^ 100 TeV 
and decays rapidly - before BBN - as expected in anomaly 
mediated supersymmetry breaking scenarios dill. Then, 
there is only constraint (ii) that the LSP's produced in grav- 
itino decay not overdose the universe, which requires y 3 / 2 < 
4x 10~ 12 (100 GcV /m LSP ). This range is indicated in Fig.[T] 
in dark gray shading, assuming ttilsp = 100 GeV. On the 
other hand, the entire range shown in Fig.^is allowed if im- 
parity is violated so that the LSP decays before BBN, if the 
LSP is much lighter than 100 GeV, if the gravitino is lighter 
than keV, or if there is no supersymmetry at all. 

It is fascinating that the correct baryon asymmetry can 
be obtained without overproducing gravitinos. The typical 
energy scale for the B violation in this scenario is about 
10 14 GeV as seen in Fig. For instance, this energy scale 
is consistent with what is expected for the Majorana mass of 
a right-handed neutrino in the seesaw mechanism [15] which 
violates B in conjunction with electroweak anomaly effects 
This would predict nearly degenerate Majorana neutri- 
nos with masses above 0. 1 eV and neutrino-less double-beta 
decay {Qv(3j3) at a rate observable in near-future experiments 
lllill . In comparison, thermal leptogenesis places a tight up- 
per bound on the neutrino mass m v < 0.11 eV Ill9ll and ob- 
servation of i)vf3j3 at larger m v would exclude that scenario. 
Another possibility for i?-violation is the D = 7 operator 
W = (UDD)(UDD)/M B , which satisfies all experimental 
constraints provided Mb <^ 100 TeV. 
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FIG. 1: The range of Trd, Mb and w that can generate an ob- 
servationally acceptable baryon asymmetry assuming a dimension 
five (n = 1) baryon number violating operator and no signifi- 
cant entropy production after decoupling is shown in dark gray and 
light gray. In supersymmetric theories, where the gravitino decays 
rapidly to an LSP with rriLSP ~ 100 GeV, the allowed region is 
restricted to the dark gray region only to avoid overclosure of the 
universe by LSP's. Both light and dark gray regions are allowed if 
the rriLSP <C 100 GeV, if the LSP decays, if the gravitino is lighter 
than keV, or if there is no supersymmetry. 



Higher Mb can be accommodated if there is significant en- 
tropy production below Trd- The regions above the curves 
in Fig. ^produce a larger baryon asymmetry and we can af- 
ford such an entropy production, which in turn also dilutes the 
gravitinos. Then, a much wider region of the parameter space 
becomes available. The predictions for a D = 5 operator, 
with entropy production below Try, are presented in Fig. [2] 

One additional concern in the case of inflation is that grav- 
itinos may be produced by quantum fluctuations in the de Sit- 
ter phase. If the gravitino mass can be ignored during the 
inflation, its coupling to the background is conformal, and no 
gravitinos are produced i^oll . Depending on the details of the 
model of inflation, the gravitino mass may be enhanced and 
hence its production |21] for helicity ±3/2 states. We have 
checked that for certain models, e.g. the supersymmetric hy- 
brid inflation model 1 22], the gravitino constraint is easy to 
satisfy [21]. Note that the helicity ±1/2 states are actually 
"eaten" inflatino which decays quickly and hence is harmless 
1E3I E4I1 . There may also be a similar concern that gravitinos 
are overproduced by brane collisions in an ekpyrotic or cyclic 
model, but this issue lies beyond the scope of this Letter. 

Finally, we point out that Af* does not have to be as high as 
the Planck scale. In fact, if the S-violation is soft, for example 
by the Majorana mass Mr of the right-handed neutrino, the 
operator in Eq. Q does not cause any unitarity violation up to 
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FIG. 2: The range of Trd, Mb and m that can generate an observa- 
tionally acceptable baryon asymmetry and avoid overclosure of the 
universe by stable LSP's expands to cover essentially both the light 
and dark gray regions in Fig. Q ;/ we allow for entropy to be pro- 
duced after the baryon-number violating interactions and gravitinos 
decouple. For each w, the allowed line in Fig.Qbecomes an allowed 
strip (shown for two cases in the Figure) because we now include the 
possibility that the baryon asymmetry may be overproduced at de- 
coupling to some degree and brought to its proper value by entropy 
production after decoupling. For the case w = 0.5, we indicate along 
the strip the values of tib /s before dilution. The additional entropy 
reduces the gravitino/LSP density, thereby opening up the allowed 
range of Trd and Mb ■ 



the Planck scale even if Ml ~ MpMp. This is an interesting 
possibility that we will not pursue further in this Letter. 

In summary, we have presented a new framework for baryo- 
genesis where CP violation lies in a gravitational interac- 
tion. The expansion of the universe promotes the microscopic 
CP violation to a dynamical violation of CPT that shifts 
the relative energies of particles and anti-particles. A CP- 
conserving _B-violating interaction in equilibrium can then 
create the asymmetry that gets frozen when the £>-violating 
interaction decouples. We have shown that it is possible 
to obtain the correct magnitude of the baryon asymmetry in 
many different cosmological scenarios: radiation-dominated 
(w w 1/3), matter-dominated (w — 0), and kinetic-energy 
dominated (1/3 < w < 1) universes. In the last case, in 
particular, one can obtain the correct baryon asymmetry while 
keeping the gravitino abundance low enough to avoid overclo- 
sure by its decay products. We envision that particle physics 
beyond the standard model can provide the required B viola- 
tion at an energy scale above 10 13 GeV, giving a whole new 
class of realistic models of baryogenesis. 
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